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ABSTRACT 

Analyses of high- redshift ultraluminous infrared (IR) galaxies traditionally use the observed optical 
to submillimeter spectral energy distribution (SED) and estimates of the dynamical mass as observa- 
tional constraints to derive the star formation rate (SFR), the stellar mass, and age of these objects. 
An important observational constraint neglected in the analysis is the mass of dust giving rise to the 
IR emission. In this paper wc add this constraint to the analysis of AzTEC-3. Adopting an upper 
limit to the mass of stars and a bolometric luminosity for this object, we construct different stellar and 
chemical evolutionary scenarios, constrained to produce the inferred dust mass and observed luminos- 
ity before the associated stellar mass exceeds the observational limit. We use the PEGASE population 
synthesis code and a chemical evolution model to follow the evolution of the galaxy's SED and its 
stellar and dust masses as a function of galactic age for seven different stellar initial mass functions 
(IMFs). We find that the model with a Top Heavy IMF provided the most plausible scenario consis- 
tent with the observational constraints. In this scenario the dust formed over a period of ~ 200 Myr, 
with a SFR of ~ 500 Mq yr~^. These values for the age and SFR in AzTEC-3 are significantly higher 
and lower, respectively, from those derived without the dust mass constraint. However, this scenario 
is not unique, and others cannot be completely ruled out because of the prevailing uncertainties in 
the age of the galaxy, its bolometric luminosity, and its stellar and dust masses. A robust result of 
our models is that all scenarios require most of the radiating dust mass to have been accreted in 
molecular clouds. Our new procedure highlights the importance of a multiwavelength approach, and 
of the use of dust evolution models in constraining the age and the star formation activity and history 
in galaxies. 

Subject headings: galaxies: evolution, high-redshift, starburst: individual: (AzTEC3) - infrared: galax- 



1. INTRODUCTION 

The origin of the dust in ultraluminous high 
redshift galaxies with redshifts z > 4 is still un- 
resolved, with supernovae (SNe), evolved AGB 
stars, the winds around active galactic nuclei 
(AGN), or growth in in o lecular clouds as possible 
sources ' ' ^ ' ' 
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origin of dust in some objects is the energy source 
powering the IR emission. An IR luminosity powered by 
luminous massive stars will suggest a high present SFR, 
favoring a SN origin for the dust. An IR luminosity that 
is powered by an active galactic nucleus (AGN) suggests 
a low present SFR and may favor the delayed injection 
of dust by AGB stars that may have formed during 
an earlier phase of much higher s tar for mation activity 
(jValiante et al.ll2009l: iDwek fc Ch crchncff l20Tll ). 

No such ambiguity exists for AzTEC-3, an ultralumi- 
nous IR galaxy located at redshift z = 5.3, t hat was first 
detected at submillimeter wavele ngths by IScott et aLI 
()2008D . It is a radio-dim source (lYounger et all I2007D 
with no X-ray counterpart (jCapak et al.ll201lD . and part 
of a young cluster, representing an overdensity of earl y 
massive galaxies at the same redshift (jCapak et al.l2011[) . 
Optical to submil limeter observations summarized in 
ICapak et al.l (|2011D show that the energy output is dom- 
inated by its far-infrared (IR) emission with luminosities 
between ~ (0.6 — 2) x 10^^ L©. The lack of any radio or 
X-ray counterparts suggests that th is luminosity is pow- 
ered by a 1800 M^m yr"^ starburst (jRiechers et al.ll2010l: 
ICapak et all [20111) . The total CO-inferred gas mass is 
~ 5 X 10^° Mq, which will be depleted in less than 
~ 30 Myr at this SFR. Such energy output suggests 
the presence of a large amount of dust, which had to 
be created (at least) in less than ~ 1.1 Gyr, the age of 
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the universe at that redshift. AzTEC-3 therefore falls 
in the category of extreme starburst galaxies with large 
dust m asses (> 10^ Mc.^) d e scribe d bv lMichalowski et alj 
(|2ninbD and lRobson etall (pOOl . 

In this paper we present a new comprehensive approach 
for studying the star formation properties and the origin 
of dust in AzTEC-3. The approach is quite general and 
can be applied to any local or high-redshift galaxy. We 
first use the UV to far-IR SED of the galaxy to deter- 
mine its bolometric luminosity and dust mass. Its SFR 
depends on the adopted IMF, and we generate a series 
of models relating the bolometric luminosity to the stel- 
lar mass and galactic age for different stellar IMFs (Sec- 
tion 2). In Section 3 wc use the observed UV-optical 
(UVO) spectrum, the UVNIR SED, and dynamical lim- 
its on the stellar mass to constrain the age of the galaxy 
from observations. Wc then generate different scenar- 
ios for the evolution of the dust, characterized by dif- 
ferent IMFs, and examine which scenario is capable of 
producing the inferred dust mass in AzTEC-3 within 
the allowed age limit, that is, before the concurrent stel- 
lar mass production exceeds the stellar mass limit (Sec- 
tion 4). Our results are discussed and summarized in 
Section 5. 

In all our calculations we adopt a flat ACDM cosmol- 
ogy, with a baryonic density parameter ilb = 0.044, a to- 
tal matter (dark-|-baryonic) density parameter of rim = 
0.27, a vacuum energy density 17a — 0.7 3, and a Hub- 
ble c onstant of Hq = 70 km s^^ Mpc~^ (jSpergel et al.l 
[200l . 

2. PROPERTIES OF AzTEC-3 

2.1. Bolometric Luminosity and Dust Mass 

Figure [T] depicts the visible to far-IR flux, vF^, from 
the galaxy calculated for a re dshift of z = 5 -3. T he ob- 
served fluxes were taken from iCapak et al.l ()2011[ ) . The 
figure shows that the energy output from the galaxy is 
dominated by the stellar radiation that is absorbed and 
reradiated by dust at IR wavelengths. Observed IR flux 
densities and their Icr uncertainties are listed in Table [TJ 
Figure [2] shows the far-IR galaxy's SED in its rest frame 
(left panel) and the observer's frame (right panel). Also 
shown in the Figure are fits of the spectra of individual 
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Figure 1. The observed visual to submm fluxes from AzTEC-3, 
calculated for a redshift of z = 5.3, are plotted as a function of 
wavelength. 



Table 1 

Summary of Submillimeter Observations of AzTEC-3^ 



Kbs (/^m) Fy(A) (mjy) a (mjy) Aq {/im)^ L^(Ao)(L0 Hz ^) 



890 8.7 1.5 141.3 1.158 

1100 7.6 1.2 174.6 1.011 

2000 3.7 1.4 317.5 0.492 

2730 0.472 ■■■ 433.3 0.0412 

3280 0.20^ ••■ 520.6 0.0173 

8165 0.087^ ••■ 1296.0 0.0116 



^ Restframe wavelengths and specific luminosities were calculated for a 
redshift of z = 5.3. 
2 3-0" upper limits. 

dust species to the observations. In order to determine 
the uncertainties in the derived dust temperatures, dust 
masses, and IR luminosities, we constructed 1000 realiza- 
tions of the observed spectrum. The flux at each wave- 
length was randomly chosen from a normal distribution 
centered around the nominal fiux. The derived quantities 
for the different dust compositions are listed in Table [2J 
The Fe and graphite dust spectra are similar to that of 
the silicate dust, and were left out of the figure for sake 
of clarity. Dust masses vary greatly from 0.3 x 10^ Mq, 
if all the dust in AzTEC-3 is in the form of amorphous 
carbon, to 3.2 x 10^ Mq, if it consists of only silicate dust. 
Because of the skin effect, the absorptivity of metallic Fe 
depends on the the grain radius, and the mass presented 
in the table was calculated for a radius of 0.33 /im for 
which the mass absorption coefficient, k, was maximum 
at the wavelength of 174.6 fim. This choice of k mini- 
mizes the mass of Fe dust required to fit the observed 
fluxes. Since AzTEC-3 is not likely to contain only pure 
Fc, silicate, or carbon-type dust, any dust consisting of 
a mixture of these elements will vary between between 
these mass estimates. 

To get a more realistic presentation of the mixture 
of dust compositions and size distribution in AzTEC-3 
we also fitted its spectrum with that of two local star- 
burst galaxies: M82, and Arp220, and a normal star 
forming galaxy NGC6946. The spectra of these galax- 
ies were obtained by using population synthesis mod- 
els to derive their stellar emission, chemical evolution 
models to derive their metallicity and dust composition, 
and radiative transfer models including the stochastic 
heatin g of dust to derive their IR spectrum and dust 
mass (jGalliano et al.l l2008f ) . Transferred to a redshift 
of z = 5.3, the IR luminosities and dust masses for all 
three galaxies were derived from their normalization fac- 
tors, obtained from the least-squares fit of their spec- 
tra to the AzTEC-3 observations. The spectra depicted 
in the Figure represent the total IR emission from dust 
residing in HII and photodissociation regions. Both 
starburst galaxies, M82 and Arp220, provided a very 
good fit to the data. Dust masses varied from 1.4 to 
2.6 X 10^ M0, and IR luminosities (~ 3 - 1000 ^im) 
from 0.79 to 1.7 x 10^'^ Lq. In contrast, the dust in 
NGC6946 is too cold, and its far-IR emission exceeds 
the 3-(T upper limits on the underlying continuum of the 
detected CO(J = 5 ^ 4) and CO(J 6 ^ 5) hnes 
(jRiechers et al.l[2010[ ) . This is evident in the inferred dust 
mass for this galaxy, which is in excess of ~ 10^° Mq. 
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Table 2 

Derived Temperatures, Masses, and IR Luminosities for Different Dust 
Compositions and Galaxy Templates 



Composition 


Td (K) 


Ma (lO»A/0) 


LiR (10"Lq) 


ft (cm^ g 


Graphite 

Silicate 

A-earbon 


32.4± 6.1 
27.5±4.5 
28.5± 4.3 
46.9±14.8 


1 4+^-2 

9 9+4.6 
— 2 2 
Q 9+2:6 
■3-^_1.4 
n Q+0.35 
^•'-'-0.17 


6+°* 

5+°-3 
U.O_Q 2 

6+°" 

U.D_Q 2 

1 l + l'^ 
^•^-0.67 


16.5 
17.6 
10.5 

28.0 


M 82 template 
Arp 220 template 
NGC 6946 template 




2.6 ±0.9 
1.4 ±0.8 

18 ± 8 


1.7±0.6 
0.8 ±0.5 
0.6 ±0.3 





^ The mass absorption coeffieient at wavelength A = 174.6 fim. 

^ Fe mass was calculated for a grain radius of 0.33 fj,m. See text for details. 
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Figure 2. The intrinsic far-IR luminosity (top panel) and the 
observed specific intensity (bottom panel) of AzTEC-3 are plot- 
ted as a function of the rest frame and observed wavelengths, re- 
spectively. Also plotted are template spectra of local galaxies and 
single-temperature fits of dust grains of different compositions to 
the observed spectrum. Data and references are given in Table [l] 

For sake of being definitive we adopt throughout this 
paper a dust mass, Mdust, of 2.0 ± 1.0 x 10^ Mq, rep- 
resenting the average value between M82 and Arp220. 
Our choice of average dust mass derived from these tem- 
plate galaxies is more realistic than any choice based on 
the assumption that all the IR emission is produced by 
a single dust specie. The average IR luminosity of these 
two galaxies is 1.0 x 10^'^ Lq and since the bolometric 
luminosity of AzTEC-3 is dominated by the IR emission 
component we adopt a bolometric luminosity, L^oi, that 
is equal to 1.0 x 10^^ Lq. 



2.2. The Star Formation Rate 

The far-IR luminosity provides a go od measure of th e 
star formation activities in galaxies (jKennicuttI [19985 • 
In a galaxy with ongoing star formation the bolometric 
luminosity and stellar mass increase with time because 
of the cumulative contribution from long-lived low-mass 
stars. Derived star formation rates depend therefore on 
the assumed star formation history of the galaxy. Fur- 
thermore, the relation between the luminosity and star 
formation rate depends on the adopted stellar IMF. Mas- 
sive stars radiate much more efficiently than low mass 
ones {L cx M^-^, in solar units) so that an IMF weighted 
more heavily towards massive stars will require a lower 
SFR for a given observed luminosity. To explore the de- 
pendence of the SFR on the stellar IMF, we used seven 
different functional forms to characterize the latter. The 
parameters characterizing the IMFs are given in Table [3l 
and Figure [3] depicts the functional form of 4 select IMFs. 
Table [3] also includes two characteristic values of the dif- 
ferent IMFs: (m), which is the IMF-averaged mass, and 
m*, which is the mass of stars that need to be formed to 
produce one Type II SN event. With these definitions, 
the stellar birthrate (yr"'^) is given by SFR/(m), and the 
SN rate by SFR/m^. 

To determine the SFR and the mass of stellar objects as 
a function of galactic age we r an the PEGASE population 
synthe sis code developed by iFioc fc Rocca-Volmerangel 
(|l997l) for the seven different stellar IMFs, listed in Ta- 
ble El for a constant SFR of 1 Mq yr"!. Table H list the 
corresponding bolometric luminosities and stellar masses 
at different galactic ages. For example, for the adopted 
Lboi = 1 X 10^^ Lq, the table shows that for a Glazebrook 
IMF the SFR is ~ 735 Mq yr'^ and the total mass of 
stars is ~ 5.6 x lO""^" Mq when the galaxy reaches an 
age of 100 Myr. At the same age, a starburst IMF will 
require a SFR of only ~ 270 Mq yr~^ to produce the 
same luminosity, and will have created ^ 7.5 x 10^ Mq 
of stars. Figure H] depicts these results in graphic form, 
showing the constant SFR required to produce a bolo- 
metric luminosity of 1 x 10^'^ Lq at a given age for 4 
different IMFs. The bold section in each line depicts the 
allowed range of galaxy ages before the model violates 
the constraint on the stellar mass, a topic discussed in 
detail in Section 3.3 below. We emphasize that this figure 
does not represent the evolution of the SFR with time. 
For each stellar IMF, it depicts the value of the con- 



Table 3 

Characteristics of the Stellar IMFs^ 



IMF 




02 


03 


Ml 


Ma 


M3 


M4 


(m>2 




Glazebrook 


1.5 


2.15 




0.1 


0.5 


100 




0.72 


70.0 


Kroupa 


1.3 


2.3 




0.1 


0.5 


100 




0.64 


89.6 


Mass Heavy 


2.35 






1.0 


100 






3.1 


52.7 


Paunchy 


1.0 


1.7 


2.6 


0.1 


0.5 


4 


100 


1.1 


66.4 


Salpctcr 


2.35 






0.1 


100 






0.35 


134.7 


Starburst 


2.35 






5.0 


100 






12.8 


24.4 


Top Heavy 


1.95 






0.1 


100 






0.78 


55.4 


^ All masses are in Mq. Stellar IMFs are given 


by: 









(f>{m) oc M-"i for Mi < M < M2; 
4i{m) oc M-"2 for M2 < M < M3; and 
<j>{m) oc M-"3 for M3 < M < M4 

The IMF averaged stellar mass. 

^ The mass of all stars born per SN event flPwek c 



CherchnefjlMTlT) 



10 



10" 



10" 



10" 




Table 4 

Starburst Luminosities and Stellar Masses-' 



Salpeter 
Glazebrook 
Top Heavy 
Starburst 



0.1 



1.0 10.0 
Mass(Mg) 



100.0 



Figure 3. The four stellar IMF discussed in detail in the paper 
as a function of stellar mass. The functional forms of these and 
the other IMFs are presented in Table [3] 



1000 



K 




100 [ 



10 



100 

age (Myr) 



1000 



Figure 4. The SFR required to produce a bolometric luminosity 
of Lbot = 1 X 10^'^ Lq at a given galactic age is depicted as a 
function of ago for four different IMFs. For each IMF, the bold 
colored line depicts the allowed range of galaxy ages before the 
model violates the constraint on the stellar mass (see Section 3.3). 

stant SFR needed to produce a bolometric luminosity of 
1 X 10^'^ Lq at the assumed galactic age. The required 



Starburst age (yr) 


10^ 




108 




109 


IQi" 


Glazebrook 














Lbol (Lq) 


9.26e4 


-09 


1.36eH 


-10 


1.74e4-10 


2.20e+10 


Mstar (Mq) 


1.03c4 


-07 


7.61e+07 


5.89e4-08 


4.50e+09 


Kroupa 














Lbol [Lq) 


7.33e4 


-09 


1.13eH 


-10 


1.54e4-10 


2.08C+10 


Mstar (Mq) 


1.05e4 


-07 


8.06eH 


-07 


6.40e4-08 


4.92e+09 


Mass Heavy 














Lbol (Lq) 


1.20e4 


-10 


1.90cH 


-10 


2.65e-fl0 


3.25C4-10 


Mstar (Mq) 


1.02e+07 


6.70eH 


-07 


3.80e4-08 


1.15e-^09 


Paunchy 














Lbol (Lq) 


7.83e4 


-09 


1.43eH 


-10 


2.05e4-10 


2.53e+10 


Mstar (Mq) 


1.05e4 


-07 


7.64eH 


-07 


4.92e4-08 


3.08e+09 


Salpeter 














Lbol [Lq) 


4.71c4 


-09 


7.49cH 


-09 


1.06e4-10 


1.49C+10 


Mstar (Mq) 


1.07c4 


-07 


8.76cH 


-07 


7.56e4-08 


6.46e+09 


Starburst 














Lbol (Lq) 


2.55e4 


-10 


3.67eH 


-10 


3.62e4-10 


3.41e-H0 


Mstar (Mq) 


9.20e4 


-06 


2.76eH 


-07 


3.17e4-07 


7.18C4-07 


Top Heavy 














Lbol (Lq) 


1.56e+10 


2.04e+10 


2.34e4-10 


2.56e+10 


Mstar (Mq) 


9.72e+06 


6.23e+07 


4.45e4-08 


3.32e+09 



^ Entries are calculated for a constant SFR of 1 Mq yr '-. 

SFR is lower for older systems because a larger fraction 
of Lhoi arises from the accumulation of lower mass stars 
over the galaxy's lifetime. 

3. THE AGE OF AzTEC-3 

3.1. Spectral Constraints on the Age 

The preceding section illustrated the importance of 
knowing the galactic age for determining the star for- 
mation, and the stellar and elemental enrichment histo- 
ries of a galaxy. The UVO spectrum can contain impor- 
tant information fo r determining this age. Figure 1 in 
ICapak et al.l (j2011l ) presents the full resolution spectrum 
of AzTEC-3, obtained with the Keck II telescope. The 
spectrum shows Si II and C IV stellar absorption features 
which are indicative of the presence of a population of 
young O and B stars. In a single burst of star formation 
the spectral contribution from young stars should dimin- 
ish with time. The presence of these features therefore 
led them to the conclusion that AzTEC-3 should be less 
than 30 Myr old. However, in an ongoing star for- 
mation scenario, the stellar population will always con- 
tain young O and B stars, and their spectral features 
will only fade at much later times as they are gradually 
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overwhelmed by the featureless continuum from the less 
massive stars. 
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Figure 5. Top panel: T he observed Keck II spectrum of 
AzTEC-3 l lCapak et al.ll2011l ) is compared to high-resolution UV 
spectra generated by Starburst99. The observed spectrum has 
been redshifted by 2 = 5.3 and rebinned to the 0.75 A resolu- 
tion of the model spectra. Model spectra were calculated assum- 
ing ongoing constant star formation rates, and are depicted for 
two stellar IMFs: Salpeter (Sp) and Top-Heavy (TH). To facilitate 
the comparison between the models and observations the spectra 
were normalized and shifted with respect to each other. The figure 
also depicts the location of important stellar lines. The observed 
AzTEC-3 spectrum is more complex since it also contains inter- 
stellar absorption lines. The main point of this panel is to illustrate 
that the observed spectrum cannot be used as an effective deter- 
minator of the age of AzTEC-3. Bottom panel: The correlation 
between the observed and the model spectrum. The colored line 
depicts the best fit line to the correlation. The synthetic spectra 
are essentially identical for all models depicted in the figure (ex- 
cept for a scaling factor), yielding an identical degree of correlation 
between the models and observations, regardless of the age of the 
system. 

To examine if the observed spectrum can be used to 
determine the age of the galaxy we constructed synthetic 
UV spectra for constant star formation scenarios and dif- 
ferent stellar IMFs (TH= Top Heavy; Sp=Salpeter) using 
Starburst99 (iLeitherer et al.lll999t iVazguez fc Leithererl 
[2001 iLeitherer et al.ll2oTo[r Figure [5] compares the syn- 
thetic spectra with the observed one. The figure shows 
that the TH model spectra have not evolved significantly 
in shape between the ages of 30 and 200 Myr. Their total 
intensity has, of course evolved, an effect not shown in 
the figure because of the offset. They are also identical 
to the Sp model at 30 Myr. The simulations show that 
when star formation is an ongoing process, the spectrum 



does not evolve significantly because even a very small 
number of 0/B stars will dominate the UV emission. 
The AzTEC-3 spectrum shows some of the same stel- 
lar absorption features present in the synthetic spectra, 
but contains in addition features caused by interstellar 
absorption, a process that is not included in the Star- 
burst99 models. The right panel of the figure shows the 
correlation between the observed spectrum and that of 
one of the models. Except for a constant offset caused 
by the evolution of the total intensity of the spectrum, 
the correlation plot is identical for all models. 

The similarity of the AzTEC-3 spectrum to the model 
spectra, and the fact that the model spectra are es- 
sentially identical even though they represent different 
epochs and were generated with different IMFs, shows 
that the observed ^ 1200 — 1600 A spectrum cannot be 
used to constrain the age of the galaxy when star forma- 
tion is an ongoing process. 

3.2. SED Contraints on the Age: The Age Degeneracy 

The UVNIR SED of a galaxy can provide independent 
constraints on its age. Using the results of the previous 
section we created a library of intrinsic stellar SEDs for 
the seven IMFs and for a grid of galaxy ages, ranging up 
to to 800 Myr. The results are shown in Figures [HI and [7] 
for 4 select IMFs and for ages between 10 and 400 Myr. 
Each panel in the figure presents the intrinsic SED (black 
line) generated by a distinct evolutionary sequence that is 
terminated at the assumed galactic age presented in the 
panel. Each evolutionary sequence is characterized by a 
constant SFR which is distinct from other sequences and 
determined by the stellar IMF and the galactic age for 
the adopted bolometric luminosity of 1 x 10^^ Lq (see 
Table g] and Figure H]). So the Salpeter 20 Myr panel 
represents the SED of a galaxy that underwent a con- 
stant SFR of - 1900 Mq yr-i for that length of time. 
In each panel, the intrinsic SED is clearly in excess of 
the observed stellar flux (red diamonds), which has been 
significantly attenuated by dust. Lacking a definitive 
dust model for the IR emission, we used the Calzetti law 
(jCalzetti et al.l[2"000( ) to characterize the extinction. For 
each evolutionary sequence we calculated the visual opti- 
cal depth, t{V), that provided the best least-squares fit 
to the observations. Values of t{V) ranged from ~ 2 — 4. 
The attenuated stellar spectra are presented by violet 
lines in the figures. 

Figure[S]depicts the minimum of the fit of the atten- 
uated spectra to the data as a function of galaxy age and 
stellar IMF. The figure shows that for each IMF there ex- 
ists an age degeneracy for which the value of is essen- 
tially unchanged. This degeneracy is caused by the fact 
that the galaxy's SED changes very little with time for 
the first ^ 100 Myr. The small changes that do occur 
are compensated for with small changes in the magni- 
tude of the attenuation. This effect is exhibited in both 
the figures of the SEDs and that depicting the value 
with galactic age. For example. Figure |8] shows that the 
Salpeter and Glazebrook IMFs provide equally good fits 
to the observed UVNIR SED for ages up to ^ 100 Myr, 
and the Top Heavy IMF for ages up to 200 Myr. The 
fits become significantly worse thereafter, an effect that 
is also clearly seen in Figures[B]and[71 The Starburst IMF 
provides an equally good fit to the observed spectrum at 
all epochs. Better fits to the observed spectra at later 
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Figure 6. Intrinsic (black curve) and attenuated (violet curve) stellar SEDs arc plotted for an adopted L^oj = 1 X 10^^ Lq, and 
for 2 different IMFs (Salpeter and Glazebrook) and 6 different galactic ages. Each panel represents a distinct evolutionary scenario, 
characterized by a constant SFR determined by the requirement that the intrinsic bolometric luminosity be equal to the adopted observed 
value of 1 X 10 Lq at the designated galactic age (see Figure|4]l. The value of Ay that provided the best fit of the attenuated spectrum 
to the observed UVNIR fluxes (red diamonds) is given in the flgurc. 
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Figure 7. Same as Figurc|6]for the Top Heavy and Starburst IMFs. 
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Figure 8. The value of for the fit of the attenuated SEDs to the 
UVNIR observations for 4 different IMFs as a function of galaetic 
age. The figure shows that for each IMF the best fitting SED is not 
sensitive to age. For ages less than ~ 100 Myr, the small changes 
in the intrinsic SED are compensated for by small changes in the 
amount of dust attenuation. The galaxy's SED provides only a 
weak limit on the age. The starburst IMF produces equally good 
fits to the observed SED for all galactic ages. 

epochs can probably be obtained if we adopted a differ- 
ent extinction law, or a time-dependent extinction law 
that reflects the changing dust composition with time. 
Such changes would need to be consistent with the far- 
IR emission as well. Such detailed investigation into the 
extinction law (i.e. optical properties of the dust, and 
the physical distribution of the dust with respect to the 
radiation sources) is currently beyond the scope of this 
paper. 

3.3. Stellar Mass Constraints of the Age: A Partial 
Lifting of the Degeneracy. 

The ambiguity in stellar age can be partially lifted by 
consideri ng the upper limits on the stellar mass in the 
galaxy. ICapak et al.l (j201lD derived a stellar mass of 
~ 10^° Mf7), usin g the Maraston library of stellar SEDs 
(jMarastonI 120051 ) and a galactic age of ~ 30 Myr. This 
age was inferred from the Keck II spectrum, assuming 
that all the star formation occurred in a single burst. 
Since we have shown that the UV spectrum cannot con- 
strain the age of AzTEC-3, we prefer to use a weaker, 
but observationally determined, constraint on the stellar 
mass. Using kinemati c CO data and l imits on the angu- 
lar size of the galaxy, IRiechers et al.l (|2010l ) derived the 
values of Mgas, the mass of the gas (dominated by the 
molecular component), and Mdyn, the dynamical mass 
of the galaxy. The dynamical mass is given by Mdyn— 
Mgas+ M stars + Mdm, where M stars and Mdm are, re- 
spectively, the mass in stars and dark matter. The dy- 
namical mass offers therefore a strict upper limit on the 
stellar mass. Adopting a gas fraction (= Mgas/ Mdyn) 
of 0.5 (jRiechers et al.|[2010D . we adopt an upper limit of 
Mstars — Mgas — 5 X IQ^^Mq ou thc stcUar mass of the 
galaxy 

Since the mass of stars builds up at different rates for 
the different IMFs, this upper limit can provide useful 
constraints on the age of thc galaxy. This effect is shown 



Figure 9. The stellar mass as a function of thc assumed duration 
of the star formation activity in AzTECS for four different IMFs. 
The dashed horizontal line corresponds to the upper limit on the 
stellar mass in the galaxy der ived from estimates o f the dynamical 
and gas masses of the galaxy l IRiechers et al.ll2010l ). For each IMF, 
the bold colored line depicts the allowed range of galaxy ages before 
the model violates the constraint on the stellar mass. 

in Figure [9l which depicts thc derived mass of stars as 
a function of the assumed age of AzTEC-3 for different 
stellar IMFs. Each symbol represents a distinct evolu- 
tionary sequence in which thc SFR proceeded at a con- 
stant rate determined by the assumed galactic age and 
stellar IMF (see Figure H]) . The Salpeter IMF requires 
the largest SFR, and is therefore the first to exceed this 
mass limit after tagei'm'O-x) 30 Myr. The Glazebrook 
IMF exceeds this mass limit after ^ 80 Myr, and the Top 
Heavy IMF after - 200 Myr. The Starburst IMF is the 
most efficient one in producing the observed bolometric 
luminosity with a low SFR, and does not have a large 
population of low mass stars. Consequently, it does not 
exceed the adopted stellar mass limit within ^ 1 Gyr, 
the approximate age of the universe at that redshift. 

4. DUST EVOLUTION: A NEW CONSTRAINT 
ON GALACTIC PROPERTIES 

The presence of (2 ± 1) x 10^ Mq of dust provides 
new additional constraints on the star formation rate 
and history of AzTEC-3. To demonstrate this ef- 
fect we used the ch e mical evolution model described in 
iDwek fc Cherchnefil ()2011f ) to follow the evolution of the 
dust for the four different astrophysical scenarios deter- 
mined by the stellar IMFs and thc upper limit of the mass 
of stars in the galaxy. For each IMF we adopted a galac- 
tic age, tageimax), taken to be equal to the last epoch be- 
fore the mass of stars produced by thc model violates the 
stellar mass limit. The ages are ^ 30, 80, and 200 Myr 
for the Salpeter, Glazebrook, and Top Heavy IMFs, re- 
spectively (see Fig. For the Starburst IMF we chose, 
quite arbitrary, a value of tage{max) = 800 Myr. The 
models then assume that star formation proceeds at a 
constant rate of - 1800, 700, 500, and 300 Mq yr^^, 
repectively (see Fig. SJ. 

We used a closed box model to follow the chemical evo- 
lution of AzTEC-3. The initial gas mass was chosen so 
that each model reproduced the inferred gas mass of 
5 X 10^° Mq at age t =tage{'max). The chemical evolution 
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Figure 10. The evolution of the dust and gas mass in AzTEC-3 are plotted as a function of time for the models eharaeterized by their 
different IMFs. Total dust masses (blue eurves) represent the sum of the dust masses produced by SNe and AGB stars. The contribution 
of AGB stars is depicted by red lines. Total and AGB dust masses are plotted for two cases: upper curves ignore the effects of grain 
destruction, whereas the lower curves adopt a Milky Way destruction efBciency. The bold section of each curve depicts the evolution over 
the time span allowed before the galaxy violates the upper limit on the stellar mass (see Fig. |9ll. The shaded region depicts the adopted 
dust mass limits inferred from the far-IR observations. 

destruction rate was calculated by integrating th e mass 
of du st destroyed as a function of shock velocity (jJonesI 
[200l . over the evolution of a supernova remnant. The 
mass of dust destroyed by a single SN remnant durin g its 
Hfetime is about m^est ~ 3 Mo pwek et all 120071 ) and 
the grain destruction rate is given by SFR x mdest/'m'*- 
The bold section of each curve depicts the evolution over 
the time span allowed before the galaxy violates the up- 
per limit on the stellar mass (see Fig. The shaded 
region depicts the range of dust masses inferred from the 
far-IR observations. 

The figure shows that when grain destruction is taken 
into account, none of the models is capable of produc- 
ing the required dust mass at any time. With no grain 
destruction, all models depicted in the figure are capa- 
ble of producing the required dust mass, albeit at differ- 
ent quantities. These idealized models can be considered 
as representatives of models in which rapid accretion in 
molecular clouds compensates for the mass of dust that 
is destroyed in the ISM. So the models show that rapid 
grain growth in clouds is an essential process for produc- 
ing the inferred dust mass in AzTEC-3. 

The model that is clearly most capable of producing 
the inferred dust mass is the one with the Top Heavy 
IMF. At the age of 200 Myr it produces more than the 
inferred upper range of the dust mass, thus allowing for 



model follows the changes in metallicities over timescales 
that are comparable to the main sequence lifetimes of the 
massive SN progenitors. Consequently, in addition to the 
delayed recycling of matter by long-lived AGB stars, the 
model also includes the delayed recycling of matter by 
massive stars. AGB yi elds for stars with ma s ses b elow 
8 Mq were taken from iKarakas fc Lattanziol (|2007l ) for 
a metaUicity of Z = 0.008. SN dust yields of massive 
stars were taken as an aver age between the zero and solar 
metaUicity yields given by iHeger fc WooslevI (|2010[ ) and 
iWooslev fc Hegeij ()2007l ). respectively (see Table[5]for se- 
lect stellar masses). Dust condensation efficiencies were 
chosen to be unity for AGB stars, and 0.5 for Type II 
SNe. 

The results of our model calculations are shown in Fig- 
ure [TOl The green curve depicts the evolution of the gas 
mass which plummets when the constant SFR depletes 
all the gas (and dust) into stars. The blue curves depict 
the evolution of the total dust mass, which is the sum 
of the contribution from SNe and AGB stars. The AGB 
contribution is a small fraction of the total dust mass, 
and is shown separately by red curves. Models were run 
for two cases: the first ignores the effect of grain destruc- 
tion, and is depicted by the solid blue curves; and the 
second takes the effect of grain destruction into account 
and is depicted by the dashed blue curves. The grain 
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Figure 11. The intrinsic UV to far-IR luminosity from AzTEC-3 
(grey line) calculated with the PEGASE stellar population code is 
plotted as a function of the rest frame wavelength, Aq. The two 
galactic templates of M82 and Arp220 represent two possible fits to 
the far-IR spectrum of AzTEC3. The intrinsic stellar luminosity of 
10-'^'' Lq falls between that of the two galactic templates (sec Ta- 
ble [2]l. The black curve represent the attenuated starlight derived 
with the Calzetti extinction curve for a value of A(V) = 2.30. 



a modest amount of grain destruction. It therefore does 
not require that the accretion in molecular clouds recon- 
stitute all the dust that was destroyed in the more diffuse 
ISM. The Starburst model is just capable of producing 
the lower range of inferred dust mass at t = 400 Myr, 
and requires that all the refractory elements returned by 
grain processing into the gas are quickly reaccreted onto 
surviving grains in molecular clouds. The models with 
a Salpeter or Glazebrook IMF fall in between these two 
models. 

Figure [TT] depicts the spectrum produced by stars 
generated by the most plausible scenario that satisfies 
both the stellar mass and the dust mass constraints of 
AzTEC-3: a galaxy with a Top Heavy IMF undergoing a 
constant SFR at a rate of 300 Mq yr~^ over a timespan 
of 200 Myr, and allowing for over half of its dust to be 
destroyed. The stellar SED is attenuated by a Calzetti 
extinction law with t{V) = 2.3. The intrinsic stellar 
spectrum is shown as the gray curve, and has a lumi- 
nosity of 1 X 10^^ Lq, which lies between the luminosity 
of the two galactic templates, M82 and Arp 220 (see Ta- 
ble [2]). This scenario is not unique, but the most plausible 
one among the models. A Top Heavy IMF was also the 
preferred scenario for the rapid and efficieii t prodtr ction 
of dust in the early universe by IGall et all (poTlal) and 
iDunne et all (|201C| ). 

5. A POPULATION III ORIGIN FOR THE DUST? 

Population HI (Pop HI) stars, the first population 
of stars being born out of a pristine, zero metallic- 
ity gas, may also be impo rtant contributors to the 
dust in these yoting galaxies (iNozawa et a l."2003l 120091 : 
iTodini fc Ferrarall2001l: iCherchneff fc Dwekii201Q). They 



are believed to be massive ( Bromm fc Larsonll2004D , and 
supported by radiation pressure at the Eddington lumi- 
nosity: 



Consider a population of massive Pop HI stars all of 
identical mass m, and let L^oi be the total bolometric 
output of this stellar population. To maintain a steady 
abundance of such stars, they have to be created at a 
rate given by: 



dN 
dt 



^bol 



^edd 



'MS 



(2) 



where tms, their main sequence lifetime, is determined 
by the time it takes to radiate all the energy released 
in the production of helium and heavy elements at the 
Eddington luminosity. It is independent of stellar mass, 
and given by: 



TMS = 0m7mc'^/L,dd « 3 X 10^ yr. 



(3) 



The SFR needed to maintain the observed luminosity is 
independent of the stellar mass and given by SFR ~ 
m{dN/dt) w 10-12 X Lboi Mq yr'^. The rate of dust 
production is given by: 



dM, 



dust 



dt 



dust 



dN 
~dt 



^dust 

tms 



Lbol 



dd 



(4) 



where Ydust is the amount of dust produced in these ob- 
jects. 

Dust formation in massive Pop HI stars has only 
been calculated for 170 M p, stars (jNozawa et al.l [20031 : 
ICherchneff fc Dwek|[20Tot ). Such a Pop HI star will have 
a luminosity of ~ 5 x 10^ Lq. The dust yield depends 
on the whether the ejecta maintains its stratified com- 
positional structure or becomes completely mixed, and 
ranges from 7 to ~ 35 Mq, for unmixed and fully 
mixed ejecta, respectively. The total mass of dust pro- 
duced at time t, assuming no grain destruction, is given 
by: 



M, 



dust 



Mrr 



(5 - 25) X — 



(5) 



where the range of ages refiects the uncertainty in the 
dust yield of such stars. The time needed to create an 
average Mdust = 2 x 10^ Mq of dust by 170 Mq Pop HI 
stars is therefore ^ 100 — 400 Myr. This timescale is sig- 
nificantly longer than their main sequence lifetime, and 
too long to maintain a pristine reservoir of gas to feed 
the formation of metal-free stars. We therefore conclude 
that the massive amount of dust seen in AzTEC-3 can- 
not have been produced by a burst of Pop HI stars that 
are currently giving rise to the observed far-IR luminos- 
ity of the galaxy. 

6. DISCUSSION 

Deriving the physical properties of high redshift galax- 
ies from their spectra is an important goal for under- 
standing their origin and evolution. In the absence of 
perfect and complete information, these derived physical 
properties are not unique, allowing for different evolu- 
tionary scenarios to explain the observations. 

In this paper we explore a new method that can, 
in principle, discriminate between different evolutionary 
scenarios. The method requires that the mass of dust 
inferred from the far-IR emission to be present in these 
galaxies be produced within an allotted time span im- 
posed by other observations. We applied our method to 
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Table 5 

Yields of Condcnsiblo Elements and Maximum Dust 
Masses in Massive Stars-' 



Element / 


20 


Mq 


50 


Mq 


Metallicity 


Z = 


Z = Zq 


Z = 


Z = Zq 


12,--; 


0.21 


0.23 


1.86 


1.61 


16 o 


1.34 


1.38 


11.2 


9.7 


24Mg 


0.065 


0.065 


0.31 


0.33 


28 Si 


0.0025 


0.19 


0.27 


0.21 


56 pg 


~ 


0.12 


0.41 


0.28 


Carbon dust 


0.21 


0.23 


1.86 


1.61 


Silicate dust2 


0.11 


0.52 


1.10 


1.00 


Iron dust 


0.0 


0.12 


0.41 


0.28 


Total dust=* 


0.32 


0.87 


3.4 


2.9 



^ Zero metallicity yields were calculated for models with 
no mixing and with explosion energies of 1.2 X 10^^ erg 
and 5 x 10^^ erg, for th e 20 and 50 Mq stars, respectively 
IIHeger & : ^yooslev^ 120101) ; Zq metallicity yields were taken 
from I WoosIeT'^nfegen II2007I ) . All entries are in units of 

Mq. 

2 The silicate dust mass was calculated assuming that all the 
Mg and Si are locked up in MgO and Si02 dust precursors. 
^ Total dust masses listed assume a 100 percent condensa- 
tion efficiency. 



the analysis of the uhraluminous submiUimeter galaxy 
AzTEC-3. The basic observations available to constrain 
the models arc the high-resolution UV-optical spectrum, 
the UV to millimeter SED, and the CO transitions that 
were detected from this object. 

A fundamental observed quantity is the bolomctric lu- 
minosity of the galaxy which in AzTEC-3 is dominated 
by the far-IR emission, and can be used to infer the SFR 
of the galaxy. However, converting L^oi to a SFR requires 
knowledge of the stellar IMF and the age of the galaxy, 
since for a constant SFR and given IMF, the bolometric 
luminosity increases with time because of the contribu- 
tion of long-lived low mass stars to the emission. 

The IR region of the spectrum can be used to derive the 
mass of dust giving rise to the emission. The dust mass 
depends on the dust composition and the dust tempera- 
ture distribution, both unknown quantities. Estimates of 
the dust mass can be derived for pure dust compositions 
(graphite, amorphous carbon, or silicates) , based on dust 
models for the Milky Way. Alternatively, the dust mass 
can also be estimated by fitting the observed SED with 
that of known galaxies for which the dust mass has been 
derived from detailed models. 

The CO emission spectrum can be used to derive the 
dynamical mass of the galaxy and the mass of molecular 
gas. Both quantities are uncertain. The dynamical mass 
estimate depends on the basic premise that the system 
has relaxed to a Keplerian disk, and the determination 
of the gas mass depends on the uncertain X-factor that 
converts the CO line intensity to an H2 mass. With- 
out any further assumptions, only an upper limit can be 
obtained for the stellar mass in this galaxy. 

The high-resolution UV spectrum can be used to deter- 
mine the age of the system. Massive young stars exhibit 
characteristic P-Cygni profiles. The presence of these 
features in the spectrum can therefore be used to deter- 
mine the age of the system. For an instantaneous burst 
of star formation, these features disappear as these mas- 
sive stars die off within ^ 30 Myr. So their presence can 
be used to estimate the age. When star formation is a 



continuous process, the massive stars are always present, 
and the strengths of the features relative to the contin- 
uum changes then much more gradually. Therefore, the 
age inferred from such spectrum will be much longer than 
that inferred for an instantaneous burst. From our sim- 
ulations, assuming a constant SFR, we find that the UV 
spectrum cannot discriminate between systems that are 
between ^ 30 to 400 Myr old, the relevant age limits for 
AzTEC-3. 

The bolometric luminosity of AzTEC-3 lies between 
0.6 to 2 X lO^'^ Lq, and to be definitive we adopted a value 
of Lboi^ 1 X lO'-'^ Lq. For the dust mass we adopted a 
value of 2 ± 1 X 10^ Mq. We adopted a very conservative 
upper limit of 5 x 10^° Mq on the stellar mass, based on 
the estimated dynamical and gas mass from this object. 

Using these basic constraints, L{,oi, M^ust, and M stars, 
we constructed models to follow the evolution of the 
galaxy's SED, stellar, and dust content as a function 
of time for 7 different stellar IMFs. Viable models were 
those that were able to produce the inferred mass of dust 
before the concurrently produced mass in stars exceeded 
the mass limit. The results of these models are depicted 
in Figure [TUl and described in Section 4. 

The results of the paper can be briefiy summarized as 
follows: 

1. Most of the dust giving rise to the far-IR emis- 
sion from AzTEC-3 must have been grown in 
molecular clouds. In the Milky Way, the need 
for grain growth in molecular clouds to explain 
the abundance of dust was first pointed out 
by [P wck & Scalo (I1980D . and more recently b y 
Zhukovska et al. (2008!) and lJones fc NuthI ()20lH) . 
Michalowski et al.. (,2010a[ ) argued that growth in 
molecular clouds is the main source of dust in high- 
redshift (5 < z < 6.5) quasars. In this paper we 
have extended this point to the ultraluminous star- 
burst AzTEC-3 using detailed chemical evolution 
models. 

2. There is a correlation between grain destruction 



12 



and the required growth in molecular clouds. At 
the peak of its dust production, which occurs at 
^ 200 Myr, the Top Heavy model can form the in- 
ferred dust mass in AzTEC-3 without resorting to 
grain growth in the ISM, but with a significantly 
reduced (compared to the Milky Way) efficiency of 
grain destruction. At earlier epochs grain growth 
in molecular clouds must account for a significantly 
larger fraction of the dust mass. In the Starburst 
model any grain destruction must be rapidly bal- 
anced by accretion in the ISM. The Salpeter and 
Galzebrook IMFs are less efficient in producing 
dust than the Top Heavy one. They can afford 
only little grain destruction, and therefore require 
a large fraction of their dust to be of interstellar 
origin. 

3. The model with the Top Heavy IMF is the most ef- 
ficient one for producing the dust in AzTEC-3, and 
therefore the more favorable model for the origin of 
the dust in this galaxy. Taking the galaxy's age to 
be 200 Myr, the longest time allowed by the stel- 
lar mass constraint, the constant SFR required to 
provide the observed bolometric luminosity and the 
required dust mass is ~ 500 Mq yr~^. The inferred 
constant SFR is significantly smaller and the age 
is significantly longer than the previous estimates 
of ~ 1800 Mq yr-^and 30 Myr, respectively, 
assuming a single burst of star formation. 

4. The Top Heavy model is, however, not unique. A 
clear discrimination between the different models 
requires a significant reduction in the uncertainties 
of the measured quantities. For example, a better 
determination of the stellar mass may have resulted 
in a lower value than the upper limit adopted in this 
paper, ruling out the models using the Salpeter or 
the Glazebrook functions to characterize the stellar 
IMF. 

5. Finally, Population HI stars could not have pro- 
duced both the currently observed luminosity and 
dust mass in AzTEC-3. 

The studies presented here illustrate the power of mul- 
tiwavelength observations and the simultaneous use of 
spectral, stellar, and dust mass constraints for determin- 
ing the properties of high-redshift galaxies. Our studies 
also highlight the need for more observational data and 
improved model input parameters. 
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